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ABSTRACT 

We have developed the first self-consistent model for the observed large-amplitude oscillations along 
filament axes that explains the restoring force and damping mechanism. We have investigated the 
oscillations of multiple threads formed in long, dipped flux tubes through the thermal nonequilibrium 
process, and found that the oscillation properties predicted by our simulations agree with the ob- 
served behavior. We then constructed a model for the large-amplitude longitudinal oscillations that 
demonstrates that the restoring force is the projected gravity in the tube where the threads oscillate. 
Although the period is independent of the tube length and the constantly growing mass, the motions 
are strongly damped by the steady accretion of mass onto the threads by thermal nonequilibrium. 
The observations and our model suggest that a nearby impulsive event drives the existing prominence 
threads along their supporting tubes, away from the heating deposition site, without destroying them. 
The subsequent oscillations occur because the displaced threads reside in magnetic concavities with 
large radii of curvature. Our model yields a powerful seismological method for constraining the coronal 
magnetic field and radius of curvature of dips. Furthermore, these results indicate that the magnetic 
structure is most consistent with the sheared- arcade model for filament channels. 



1. INTRODUCTION 

High-cadence Ha observations of large-amplitude 
standing o scillations a l ong a filament axis were first re- 
ported by Uing et ail (|2003[ ). who measured bulk dis- 
placements of up to 140 Mm, with a period of 80 min 
and a high velocity amplitude of 92 km s -1 . Promi- 
nence oscillations are classified as "large amplitude" 
when velocity is of the order of 20 km s _1 or larger 
(jOliver fc Ballester 2002). Similar large-amplitude lon- 
gitudinal ( LAL) oscill ations were observe d later by 
IJing et al.l ([2006h and IVrsnak et all ((20071 ). All the 
LAL oscillations reported so far show periods of 0.83 — 
2.66 hours, damping times of 2 — 10 hours, damping time 
per period of 2.3— 6.2, and amplitudes of 30— 100 km s . 
These oscillations apparently were triggered by an ener- 
getic event: a sub- flare, a microflare, or a flare close to 
the filament. In addition, prominences undergo large- 
amplitude transverse oscillations in response to waves 
emanatin g from dis tant flares (e.g., Gilbert et al.l [20081 : 
Tripathi et al. 2009). The work presented here addresses 
only the longitudinal motions. 

Several m odels have been p roposed to explain the LAL 
oscillations. IJing et aT] (|2006 ) interpreted their observed 
oscillations in terms of the IKleczek fc Kuperusl (|1969f ) 
model. In this model the prominence is considered a 
slab with a magnetic field along the filament, the restor- 
ing force is the magnetic tension, and the motion is per- 
pendicular to the main axis of the filament. However, 
this geometry is inconsistent with the observed direc- 
tion of the LAL motion and the most likely magnetic 
structure: dipped or twisted field nearly parallel to the 
polarity inversion line (PIL) (Mackay et al. 2010). The 
observations show that the displacements of the filamen- 
tary structures comprising the prominences are mostly 
aligned with the PIL and hence along the filament axes, 
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without visible changes of the magnetic structure af- 
ter the onset of the oscillation. In contrast, the mag- 
netic tension drives only transverse motions , which would 
be perpendicular to the LAL oscillations. IVrsnak et al.l 
(2007) assumed a flux-rope geometry, and suggested 
that these oscillations are driven as a longitudinal-mode 
standing wave on a tightly coiled spring with fixed ends. 
This model also predicts motions perpendicular to the 
local magnetic field, inconsistent with the observed LAL 
dynamics. Another possibility is that the restoring force 
is associated with gas pressure differences. However, the 
temperature differences needed to accelerate the threads 
to the observed speeds are on the order of several mi l- 
lion Kelvins, which are not observed ((Vrsnak et al.l2007h . 
The damping mechanism is also poorly understood. Sev- 
eral damping mechanisms have been suggested but not 
rigorously tes ted, e.g., energy leakage b y the emission of 
sound waves (Klcczck & Kupcrus 1969) or some form of 
dissipation (|Tripathi et al.ll2009l : IOliverll2009D . 

In this work we propose explanations for the forcing 
and damping of the LAL oscillations based on our pre- 
vious inve stigations of prominen ce formation and evolu- 
tion (e.g.. lAntiochos et al.l [20001. We found that all of 
the resulting condensations oscillate after their forma- 
tion, and that these motions damp after several oscilla- 
tions. By studying in detail the motion and damping 
of many simulated thread s in a sheared-arcade filament 
channel (|Luna et al.ll2012l ) , we could identify and analyze 
quantitatively the underlying physical processes. Our 
prominence model thus provides the first self-consistent 
explanation for the LAL oscillations. 

2. NUMERICAL MODEL AND RESULTS 

We focus our attention on the 47 individual oscil- 
lati ng threads model ed in our recent work (for details 
see iLuna et ail 12012ft , which form withi n dipped flux 
tubes of a sheared-arcade filament channel ([DeVore et al.l 
2005); all flux tubes are assumed to remain rigid. The 
one-dimensional equations governing the plasma dy- 
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namics and energetics in these tubes are numericall y 
solved with our ARGOS code (jAntiochos et al.lll999D . 
In response to steady, localized heating imposed in the 
tube footpoints, condensations of cool, dense plasma 
are produced in the corona b y the well-investi gated 
thermal noncquilibrium process (Antiochos & Klimchuk 
1991 1 : JAntiochos et all [1 999. 20 001: IKarpen et al.1 l200lT 
2003 , [2005, 20061: IKarpen fc AntiochosT BroOS: Xi alataLl 
2011 . 

We found that the combination of heating asymmetry 
and tube geometry determines the (usually off-center) 
position where the condensation forms, which is always 
closer to the footpoint with weaker heating, and the 
initial velocity of the thread. Threads rapidly settle 
down after oscill ating around the bottom of the dip 
(jAntiochos et all 12000 : IKarpen et all 12001 l200l . and 
accrete mass thereafter at a constant rate due to steady 
footpoint heating. As discussed in detail below, we find 
that the oscillatory features of the model prominence 
threads are strikingly similar to those of the LAL os- 
cillations. The threads move mainly along the filament 
axis because the local magnetic field is highly sheared 
(i.e., nearly aligned with the PIL). 
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Fig. 1. — Plot of the temporal evolution of the bulk velocity of 
the center of mass of the condensation example (solid line). The 
decaying exponential fit for the strong damping (dashed line) and 
weak damping (dot-dashed line) are also plotted (see text). Ad- 
ditionally, the theoretical velocity fit obtained with our pendulum 
model described by Eq. \5\ is shown (dotted line) . 

A typical example of an oscillating thread is shown in 
Figure [TJ Around t c = 2.7 hours, the thread forms with 
an initial velocity amplitude of 49 km s _1 but rapidly 
reaches a maximum speed of 53 km s . The maximum 
displacement in this case is 40 Mm with respect to the 
bottom of the dip (the approximate equilibrium posi- 
tion). We note that the oscillation is very regular in 
time with a constant period of 1.0 hour. The oscilla- 
tion is damped, but the velocity profile cannot be fitted 
with a sinusoidal function with a single decaying expo- 
nential envelope throughout the entire simulation. The 
initial oscillations damp rapidly at first followed by much 
weaker attenuation. We have fitted an exponential func- 
tion for the first three periods of the oscillation enve- 
lope (dashed line) , yielding strong damping with a decay 
time of r s = 2.3 hours. We have also fitted a decaying 



exponential function for the envelope of the curve after 
t = 10 hours (dot-dashed line), yielding a weak damp- 
ing time of r w = 13.1 hours. The existence of these two 
damping times indicates the existence of two physically 
distinct phases in the temporal evolution of the thread. 
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Fig. 2. — Scatter plot of the velocity amplitude as a function 
of the displacement amplitude for the set of 47 oscillati ng threads 

i filled circles), a nd observational data (diamonds) from Ji ng et all 
2001 . |20M) and IVrsnak et ahl J2007D . 

Similarly, we have computed the oscillation parame- 
ters for the 47 threads identified in our simulations. Fig- 
ure [5] indicates that a higher displacement amplitude is 
associated with a higher thread velocity. Most of the 
thread oscillations have velocity amplitudes larger than 
20 km s _1 , and hence are classified as "large amplitude" 
oscillations. We have not found any dependence of the 
periods on the length of the flux tubes or the mass of the 
condensations. However, we have found a dependence on 
the average radius of curvature of the dip, R, as shown 
in Figure [3] In <j3] we explain the reason for such a de- 
pendence. There are few cases that do not fit very well 
the model curve. These tubes have a irregular shape and 
the radius of curvature is not uniform in the dip. The 
strong damping time r s as a function of the weak damp- 
ing time r w is plotted in Figure @] Clearly t s < t w even 
in the cases where the weak damping time is small. The 
decay time per period is t s /P =1—6 for the strong and 
t w /P = 2-27 for the weak. 

The clear agreement between the observed range of 
amplitudes, periods, and damping times ( JQ) and the re- 
sults of our simulations (see Fig. [2]) led us to speculate 
that the reported LAL oscillations could be described 
with our theoretical model. We propose the following 
scenario to explain the observed periodic displacements, 
quasi-aligned with the filament axis, of parts of the fila- 
ment. The LAL oscillations are triggered by a n energetic 
pertu rbation at one end of the filament (e.g.. lJing et al.l 
120061 ): the nature of the trigger is discussed in S}5] This 
perturbation pushes the threads to oscillate at the same 
time, so the LAL oscillations represent the collective mo- 
tion of bundles of threads that initially move in phase. In 
our model, the threads oscillate along the dipped part of 
magnetic flux tubes that are more or less aligned with the 
axis of the filament. Only a few LAL cycles were observed 
after the trigger, which we associate with the strong ini- 
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tial damping of the oscillations of our simulated threads. 
The weak damping has an small contribution to the at- 
tenuation in the hrst cycles. In addition, after several 
oscillations the excited threads lose their coherence due 
to slight differences in the oscillation periods, explaining 
why different parts of the fi lament oscilla t e incr easingly 
out of phase as observed by iVrsnak et al.l (|2007t ) . 
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Fig. 3. — Scatter plot of the oscillation period of the threads as a 
function of the mean radius of curvature of the dipped part of the 
tubes (filled circles). The solid line shows the theoretical period 
from the pendulum model described by Eq. [4] 
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Fig. 4. — Scatter plot of the strong damping time r s as a function 
of the weak damping time r w ■ 



3. OSCILLATION PHYSICS: RESTORING FORCE AND 
DAMPING 

We noted the existence of thr ead oscillations in pre- 
vious numerical s t udies (e.g., lAntiochos et al.l l2000t 
iKarpenet all 120011 I2003D . but did not explore the ori- 
gin of the restoring force and damping. Two possible 
origins for the restoring force have been proposed: pres- 
sure gradients and gravity. For several reasons, we find 
that pressure forces cannot account for the observed or 
simulated behavior. Small-amplitude prominence oscil- 
lations were studied in two-dimensional mag netic config- 
urations with slight flux-tube curvatures by I Oliver et al.l 



(1992 ). IQliver et al.l (fl99l. IJoarder fc Roberts! (fl99l . 
and lOliver fe Ballesterl (|1995l) . who concluded that the 
oscillations were pressure driven. However it is un- 
clear whether these linear studies are applicable to large- 
amplitude, low-frequency oscillations. In the pressure- 
driven case, rapid changes in threa d mass and length 
would cause frequency shifts (e.g., iMorton fc Erdelvll 
l2009t lErdelvi et al.ll20lH ) that are not observed in either 
the solar or the simulated LAL oscillations. Finally, our 
simulations reveal that the pressure force on the con- 
densation is negligible in comparison with the gravity 
force. Therefore, we conclude that pressure-driven oscil- 
lations are inconsistent with both the LAL observations 
and with the behavior of our simulated threads. 

We propose instead that the restoring force is mainly 
the solar gravity, and the damping mechanism is mainly 
the ongoing accretion of mass by the condensation. In 
the dipped part of a flux tube, the projected solar grav- 
ity is always directed towards the bottom of the dip. 
Therefore the thread oscillation resembles the motion of 
pendulum with a length equivalent to the average radius 
of curvature of the dip, R. The temporal evo lution of 
the cool mass (see Fig. 3a of iLuna et al.ll2012h consists 
of an initial transient phase of strong growth starting at 
time t c and a linear regime after to. The transient phase 
is short (to — t c < P/2) so the mass of the condensation 
is well described as 



m(t) = too + a (t — to) , 



(1) 



where a is the growth rate of the condensation mass, 
and too is the mass at the beginning of the linear regime 
at to. In addition, all parts of the condensation move 
with similar velocities and the motion is well described 
by the center-of-mass velocity. Although the threads pe- 
riodically expand and contract slightly throughout the 
simulations, these perturbations do not affect the motion 
of the center of mass. We also assume that the ampli- 
tudes of the thread displacements are small with respect 
to the curvature radii of the dips, consistent with our 
simulations. The condensation exchanges mass and mo- 
mentum with the surrounding corona, thus introducing 
the effects of drag into the force balance. Thus we obtain 
the equation of motion as 



and 



d 2 S 

dt 5 " 


+ t 


1 ds 
-t' dt 


+ u> 2 


s = 


t' = 


= *0 


too 
a 


= *o- 


"00 

u 



0. 



(2) 



(3) 



where s(t) is the position of the center of mass of the 
condensation along the tube, lo is the angular frequency 
of the oscillation, and the phase ipo = u> mo/a. Equation 
^ is similar to the equation of a simple gravity-driven 
pendulum, with an extra drag term (second term). The 
coefficient of this drag term is the mass growth rate per 
unit of mass, ,.■, . For the initial stages of the evolution 

' m{t) & 

(t ^ £') the drag is intense, but for times much larger 
than t' this term is negligible. 
The angular frequency of the oscillation is 
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(4) 
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where P is the period and go is the gravitational acceler- 
ation at the solar surface. Equation ([2j is the zero-order 
Bessel equation with solution Jo[uj{t — t ) + tpo]. This 
oscillating function describes very well the first stages of 
the evolution: the amplitude decreases with time but not 
exponentially, which explains why we cannot fit a sim- 
ple decaying exponential to the velocity profile shown in 
Sj2] After several periods the Bessel function is basically 
a sinusoid. However, the simulations also exhibit much 
slower damping during the later stages of the oscillation 
(see Fig. [T]) , so it is necessary to introduce an additional 
damping term. The thread displacement takes the form 



equivalently 



s(t) = s Q + A J [oj (t - t ) + ^ ]e- ( *- to)/n " 



(5) 



The velocity derived from this equation is over-plotted 
in Figure []] (dotted line). We obtain a very good 
match between the simulated oscillation and the oscil- 
lation model, which is also valid for the oscillations of 
the other 46 simulated threads. The additional weak 
damping might be associated with non-adiabatic effects 
such as conductio n and radiation (jTerradas et al.l l2001t 
ISoler et all |2009() . These losses are very effective in 
straight flux tubes where the curvature and gravity are 
not included, and the restoring force is the pressure gra- 
dient. In those systems, the non-adiabatic losses weaken 
the restoring force in each oscillation, damping the move- 
ment of the thread. However, in our system the main 
restoring force is the gravity projected along the tube, 
and the non-adiabatic losses are not so effective. 

Equation (JS|) indicates that both mass accretion and 
non-adiabatic losses contribute to the strong damping at 
the beginning of the oscillation. However, the former 
is more important because the latter term has a long 
damping time. The effective decay time associated with 
the mass loading is implicit in the Bessel function because 
its amplitude decreases rapidly in the first few cycles. 
This damping is related to the phase ipo of Equation ([5]) . 
For i/}q = the damping is maximum and the damping 
time takes its shortest value, but for larger values of ipo 
the damping is weak and the Bessel function is basically a 
sinusoid. Because this phase is inversely proportional to 
the mass accretion rate a (see Eq. [3]) , rapid mass loading 
implies strong damping of the oscillation and vice versa. 
However, relatively large values of the loading rate imply 
that tf>o p» 0, indicating that the effective damping time is 
not very sensitive to a variations. The damping depends 
also on mo , which is not straightforward to characterize 
with a simple analytic approach. For these reasons it is 
difficult to find a simple relation between the damping 
time and the parameters of the tube and the heating. 
This relation will be the subject of a future study. 

4. PROMINENCE SEISMOLOGY 

In our prominence model the magnetic structure exists 
independent of the condensations. The shear and bend- 
ing of the coronal magnetic field are determined by the 
shear at the photosphe re and the tension of the o verly- 
ing field along the PIL (|DeVore fc Antioch os 2000). The 
magnetic-field structure is self supporting, and the in- 
teraction with the plasma content is small. However, 
the magnetic structure must support the denser threads. 
Thus, the magnetic tension in the dipped part of the 
tubes must be larger than the weight of the threads or 
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where B is the magnetic field strength at the bottom of 
the dip, n is the particle number de nsity in the thread , 
the mean particle mass m = 1.27 m p (jAschwandenl2 004) . 
and TOp is the proton mass. Combining Equations (jj) 
and d6| we obtain 
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which constrains the minimum field strength as a func- 
tion of the thread oscillation period. Equation ([7]) can 
be written as 



BIG] > 26 



1/2 



P [hours] 



(8) 



For the thread electron number densities of our simula- 
tions (n = 1.6 — 6.2 x 10 11 cm~ 3 ) and the oscillation pe- 
riods (P = 0.75 — 1.75 hours), Equation ((HJ yields a min- 
imum field strength of B > 31 — 75 G for the 47 threads 
with a mean value of 51 G, in agreement with the sparse 
observ ations of intermedia te-type prominence magnetic 
fields (Mackav et al.l 1201 01 ). In our prominence model 
the deep-dip flux tubes have an average area-expansion 
factor of 3 with respect to their footpoints. Thus we 
estimate the mean field strength at the flux-tube foot- 
points as 153 G, a lso consistent with observations (e.g., 
IAschwande"n1l2004D . 

In addition to providing valuable information on the 
difficult-to-measure coronal magnetic field in filament 
channels, our model also can determine the radius of 
curvature of the prominence flux-tube dips from the 
observed oscillation p roperties. We ha ve used the 
Jing et al. (2003) and Vrsnak et al. (2007) observations 
to solve Equation @ for R. The observed periods are 
80 min and 50 min respectively, yielding corresponding 
curvature radii of 152 Mm and 62 Mm. This large-scale 
curvature of the field deduced from the observed oscilla- 
tions, combined with the large observed thread displace- 
ments, argue strongly for the applicability of the sheared- 
arcade magnetic structure. Assuming that the density of 
the observed prominence threads is in the range of the 
values found in our simulations, we estimate the mini- 
mum field strength in the prominence as approximately 
34 — 68 G and 23 — 46 G, respectively, for the two cases. 

5. THE OSCILLATION TRIGGER 

In our simulations, steady footpoint heating produces 
condensations that oscillate after formation. Each con- 
densation has an initial velocity imparted by the ther- 
mal nonequilibrium process and determined by the heat- 
ing asymmetry and flux tube geometry. However, the 
observed LAL oscillations clearly are related to an ener- 
getic event (subflarc, microflare, or flare) near the fila- 
ment ends. Although a detailed study of the excitation 
mechanism is beyond the scope of this letter, we propose 
that the energetic event increases the heating impulsively 
at one end of the filament, increasing the evaporation of 
chromospheric plasma into the affected flux tubes. Our 
preliminary studies have revealed that condensations are 
remarkably robust; increased heating at one footpoint 
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produces a flow of hot evaporated plasma that pushes 
a thread to a new location, rather than destr oying it in 
situ. For example, in the Uing et all (|2003f) event the 
subflare pushed an entire section of the filament 70 Mm 
away from the initial position. We speculate that extra 
heat is deposited only on those footpoints located clos- 
est to the energetic event, so only parts of the filament 
are excited simultaneously and oscillate collectively in 
phase. In addition, the increased evaporation produces 
more mass accretion and increases the oscillation damp- 
ing. 

6. CONCLUSIONS 

We conclude that the observed LAL motions represent 
a collective oscillation of a large number of cool, dense 
threads moving along the magnetic field, triggered by 
a nearby energetic event (subflare, microflare, or flare). 
The main restoring force is the projected gravity in the 
flux tube dips where the threads oscillate, even in threads 
with large radii of curvature. The period is independent 
of the tube length and of the constantly growing mass. 
The oscillations are strongly damped by the mass accre- 
tion of the threads. After several oscillations the excited 
threads lose their coherence due to slight differences in 
the oscillation periods, causing different parts of the fil- 
ament to oscillate increasingly out of phase as observed. 



The mechanism that initiates the collective oscillation is 
not determined in the present investigation. We spec- 
ulate that the energetic event impulsively increases the 
evaporation at the filament flux-tube footpoints closest 
to the energetic event, increasing the evaporation rate 
only on one side of the threads. The hot evaporated 
plasma pushes the existing threads, initiating the oscil- 
lations. A detailed study of the initiation mechanism 
will be the subject of future work. More observations of 
oscillating prominences are needed to further refine and 
test our model. 
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